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I. Abstract 

In earlier investigations, the adaptation and implementation of a modified two-level cor- 
rections (or targeting) process as the onboard targeting algorithm for the Trans-Earth In- 
jection phase of Orion is presented. The objective of that targeting algorithm is to generate 
the times of ignition and magnitudes of the required maneuvers such that the desired state 
at entry interface is achieved. In an actual onboard flight software implementation, these 
times of ignition and maneuvers are relayed onto Flight Control for command and execution. 
Although this process works well when the burn durations or burn arcs are small, this might 
not be the case during a contingency situation when lower thrust engines are employed to 
perform the maneuvers. Therefore, a new model for the two-level corrections process is for- 
mulated here to accommodate finite burn arcs. This paper presents the development and 
formulation of the finite burn two-level corrector, used as an onboard targeting algorithm 
for the Trans-Earth Injection phase of Orion. A performance comparison between the im- 
pulsive and finite burn models is also presented. The present formulation ensures all entry 
constraints are met, without violating the available fuel budget, while allowing for low-thrust 
scenarios with long burn durations. 


II. Introduction 


NDER nominal operational conditions, the Crew Module’s (CM) 33,361 N main engine 
U performs the 3-burn Trans-Earth Injection (TEI) sequence. However, the spacecraft must 
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also have the ability to autonomously target and execute the maneuvers to return the crew 
safely to Earth using the backup auxiliary engines. Since the total thrust of these eight (8) 
auxiliary engines is 4,448 N, the burn durations are naturally much longer. For instance, the 
baseline TEI-1 maneuver, implemented by the main engine, lasts roughly 5.5 minutes. In 
contrast, if TEI-1 is executed using solely the auxiliary engines, the duration increases to 55 
minutes. Due to these extended burn durations, it is no longer accurate to approximate each 
maneuver as impulsive. Thus, it is necessary to develop an autonomous targeting algorithm 
that allows maneuvers of finite duration during a main engine failure scenario. 

Autonomy, for the Orion trans-Earth phase, refers to the ability to (a) automatically 
identify a suitable startup arc 1-3 and (b) use that solution to successfully target the specified 
entry constraints within the fuel budget available at the time. The first step, the identifi- 
cation of the startup arc, can be accomplished in one of two ways. The simplest and most 
common approach is to generate a database of optimal solutions over a time interval of in- 
terest and use those as “nominal” departure scenarios at the desired time. 1,2 The targeting 
process then reconverges the solution as needed to account for discrepancices in the timing 
and state. More recent methods 3 consider the use of infeasible solutions (i.e. with state 
and time discontinuities) based on a series of two-body approximations. Both methods are 
suitable for the generation of an initial guess in this case. However, from a historical per- 
spective, the database method has been successfully employed since the Apollo era, though 
more commonly from a ground operations perspective. In an onboard determination sce- 
nario, the database method allows for reduced computation time when the database includes 
sample optimal solutions at an adequate rate. Problems that are time sensitive, such as the 
Orion TEI sequence, require an increased number of samples, roughly one every 12 hours. 
The examples presented here employ the database approach to extract an initial guess for 
the subsequent targeting process. The initial guess supplied to the targeting process does 
not meet the specified path constraints, and sometimes the solution may not meet the cost 
constraint (i.e. fuel available). The solutions supplied are also based on the availability of 
the main engine. As such, in a main engine failure scenario, the quality of the initial guess 
supplied degrades significantly. The present study is strictly focused on the second stage 
of the autonomous targeting process, re-targeting the entry interface state using only the 
resources available onboard at the time (i.e. fuel left and operational engine) based on the 
initial guess supplied. 

Significant research has been done on the subject of optimal finite thrust guidance. 5-8 
Among these methods, nonlinear programming is commonly employed in solving optimal 
and nonlinear targeting problems. 7,8 The process of identifying, numerically, optimal or 
feasible solutions via nonlinear programming is basically the same. The main difference 
is that optimization problems require a cost index be specified and feasibility problems, 
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such as constrained nonlinear targeting, do not. Of course, the identification of feasible 
solutions that meet all the specified constraints is also accomplished through linear targeting 
methods . 4,9 These classical methods employ the state transition matrix to compute the 
necessary constraint gradients during the corrections process. More recent studies 10 employ 
a similar approach to compute “analytic” derivatives for implementation in a nonlinear 
programming process for trajectory optimization. Naturally, a nonlinear process is preferred 
when the computational resources are available. However, for onboard determination, the 
optimality of a solution is not as critical as the availability of a feasible solution. In this 
case, the inherent simplicity of linear targeting algorithms leads to a reduced cost in flight 
software development and validation. 

Earlier studies consider optimization methods for use during onboard targeting processes. 
These include the use of a simplified adaptive guidance law for targeting relative to a 
predetermined nominal trajectory 11 or implementation of an efficient sequential gradient- 
restoration algorithm employing multiple subarcs . 12 These studies, though, are tested for 
orbital transfer and rendezvous, which do not have the third body effects that so greatly 
impact the mission in this study. The algorithm presented here is a modified two-level cor- 
rector , 4,9, 13-16 employed during the design of the Genesis trajectory , 9,14,15 that allows for the 
incorporation of finite burn maneuvers. 

A two-level targeter (or corrector) is primarily based on linear system theory; it uses a 
time-varying linearized dynamical model and a minimum norm solution to compute solu- 
tion updates. These linear updates are implemented in the nonlinear system in an iterative 
corrections process that repeats until a feasible solution is identified in the vicinity of the 
startup arc. The two-level process offers several advantages: because the updates are based 
on the linearized model, it is numerically simple and computationally efficient. It does not 
require knowledge of a nominal solution, relying instead solely on the current path of the 
vehicle. The two-level correction process also allows for straightforward addition of path 
constraints, both those at specific points (i.e. entry interface ) 4,9,15,16 and those applied over 
the trajectory as a whole . 16 However, it was originally designed to use impulsive maneu- 
vers as control variables. In this investigation, the classical impulsive two-level corrections 
process 4 is modified to incorporate accurate thruster models to allow for burns of hnite du- 
ration while still retaining the structure and simplicity of the original algorithm so that it is 
suitable for onboard calculations. The theoretical elements of the formulation are presented 
next, followed by a series of performance comparisons between the impulsive and hnite burn 
targeters. For all the comparisons, the same startup arc (based on impulsive maneuvers) 
and the same entry targets are used. 
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III. Finite Burn Targeting Algorithm 


The basic structure of the Level II process in a two-level targeter that incorporates finite 
burns is the same as that of the impulsive two-level targeter. 4 First, both algorithms treat 
the position and times of each patch point as control variables. Furthermore, both algorithms 
employ a minimum norm solution in computing the updates to these control variables. The 
differences, which subsequently lead to added complexity and computational overhead, stem 
from the increased dimensionality of the state vector associated with any burn arc. Due 
to the interdependency between these state variables, the partial derivatives are also more 
complex in nature than those of the impulsive targeter. 

Traditionally, an impulsive two-level targeter requires a startup arc represented by a 
series of “patch states.” These states, also termed “patch points,” are selected by the user 
as representative waypoints along the trajectory. The user supplies the time and state at 
each patch point, t k and vjl" ] for k — 1, ■ ■ ■ , IV, respectively. Each state xjt is 

then numerically integrated forward over an interval [t kl tk+ 1 ], for k — 1, • • • , N — 1. The 
integrated state, at time t k + 1 , is recorded as xT_ 1 . This is to allow for the possibility that 
the user supplied velocity at that point, v^ +1 , may not coincide with that identified during 
the propagation, vjT +1 . Such differences may arise due to a previously scheduled impulsive 
maneuver at that point or to differences in the models used (two- vs. three-body). This is 
graphically illustrated in Figure l(a)-l(b). Thus, a Level I process leads to a trajectory that 
is continuous in position but, potentially, discontinuous in velocity at certain points. This is 
rectified by incorporating a Level II correction. 

The Level II process adjusts the positions and times of each free patch state to drive any 
of the interior velocity discontinuities to zero, as well as meet any additional user specified 
constraints. This is graphically illustrated in Figure 2(a)-2(c). Figure 2(a) is representative 
of the scenario in Figure 1(b). Figure 2(b) illustrates how the patch state positions, and 
potentially the associated times, have been adjusted by the Level II process. Subsequently, 
since the corrections are linear in nature, 4 propagation of the updated patch states in the 
nonlinear system can lead to a trajectory that is, once again, discontinuous in position. The 
Level I process is sequentially applied once more to generate an updated trajectory that is 
continuous in position. The combined Level I and Level II processes are generally repeated 
until the user specified tolerances are met for position and velocity continuity, as well as any 
additional constraints specified. Additional constraints may include velocity continuity at 
all patch states except where maneuvers are allowed, and interior or boundary constraints, 
among others. 4 

It is important to note that the initial guess need not be feasible. That is, posi- 
tion/velocity/time continuity is not necessarily required for the targeter to successfully con- 
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(a) Impulsive Two-Level Targeter: Before Level I Process is Applied 



(b) Impulsive Two-Level Targeter: After Level I Process is Applied 

Figure 1. Level I Process 


5 of 28 


^fc+l> r k + 1 



tfc-h r k — 1 

(a) Impulsive Two-Level Targeter: Before Level II Process is Applied 


r ' 

L k+V ' k+1 


t'. 


k—V k—1 


t' J 
r k 


(b) Impulsive Two-Level Targeter: Corrections Suggested by Level II Minimum Norm Solution 



(c) Impulsive Two-Level Targeter: After Propagating the Level II Updated Patch States 


Figure 2. Level II Process 
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verge. However, since the overall process is based on linear systems theory, the initial 
discontinuities can impact the computation time. An initial guess with large discontinuities 
leads to an increased number of iterations. Naturally, an initial guess with absurdly large 
discontinuities can lead to non-convergence. Of course, a low quality initial guess can have a 
negative impact on both linear and nonlinear targeting algorithms. However, linear targeters 
will naturally be more sensitive to large errors. Developing a good initial guess is a problem 
within itself and highly dependent on the particular application of interest. 

Provided a suitable initial guess is available, the formulation of the impulsive two-level 
targeter 4 is generalized in nature. As such, it can be applied to any problem that employs 
impulsive corrections. However, problems that employ continuous control of any kind cannot 
benefit from this approach, at least not in its original form. The key to transitioning the 
methodology to address problems that include segments of continuous control is to formulate 
the control variables in terms of constant parameters that can be adjusted. For example, if 
the thrust vector is inertially fixed, and the engine only allows fixed thrust or acceleration 
levels, the control variables become the time of ignition, and the direction and duration of 
the burn. Under similar conditions, if linear steering is allowed, the control variables become 
the time of ignition, the duration of the burn, the initial burn direction, and the rate of 
change of the burn direction. 

In the classical impulsive two-level targeter, the Level I process employed Av’s at the 
start of each segment to achieve position continuity. These Av’s, and - if desired - the time 
at which the maneuvers are executed, are control variables in that case. In a Finite burn 
process, the Level 1 control variables include the ignition time, burn time, and thrust vector 
parameters. The structure of the finite burn two-level targeter is subsequently developed 
and presented here. 

A. Level I Process 

As previously discussed, the application of a Level I process 4 to the orbital transfer problem 
typically involves the identification of an arc that spatially connects two points in space. This 
is the n-body equivalent of a two-body Lambert targeter, except the time of flight is not 
necessarily fixed or pre-specihed. This ultimately reduces to some form of linear differential 
correction where Av’s are adjusted to meet the specified goals. In the present study, however, 
impulsive maneuvers do not adequately model the true nature of the burn implementation. 
Thus, the Level 1 process traditionally employed in the two-level targeter 4 requires some 
modification to incorporate finite burn arcs. 

Consider a segment defined by patch points k — 1 and fc, as shown in Figure 3. In a 
Level 1 process that employs finite burns rather than impulsive maneuvers, the burn arc is 
considered to be a subsegment of the arc between patch points k — 1 and k. The end of 
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the burn occurs at point T. In identifying finite burn arcs, it is necessary to consider an 



Figure 3. Level 1 Process 
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le thrust direction associated with patch 


r k v fc m k m 9k u k 
spacecraft mass, the propellant flow rate, and t 
point k, respectively. The goal is to identify a relation between the target, which is the 
terminal position vector at point k (r fc ), and the control variables. The control variables are 
the thrust direction (ufc_i = u*,) and the time at the end of the burn (tx)- The variational 
equation for the burn subsegment is, 


5r T — v T dtT 


<5r fc _i - v^Stk-i 

Svrp a T StT 


Sv t-i - a+^Stk-i 

5m ? + rh~ T 5t T 

= $(T,fc- 1) 

5ml_ x + rh+ k _i6t k - 1 

5m~ — m~ 5tx 


5m+ k _i - m+^Stk-i 

hUy — U 7p5tT 


Su t-i ~ 


for the burn subsegment, where $(T, k — 1) is the state transition matrix between patch point 
k and point T. As in the impulsive formulation, the state transition matrix is partitioned 
into sub-matrices corresponding to each state: 


$(T,fc- 1) 


A-T,k- 1 

Bx,k - 1 

pT,k-l 

Fr,k - i 

Gr,k - i 

Cr,k - i 

Dr,k - i 

Hr,k - i 

It, k-i 

Jr,k-i 

Kr,k-i 

Lr,k - i 

Mr,k - i 

Ar,fe-i 

Or,k-i 

Pr,k - i 

Qr,k - i 

Rr,k - i 

Sr.k-\ 

Tx,k - i 

Utm-\ 

Vr.fc - i 

Wr,k - i 

Xr,k - i 

Yr,k-i 


( 2 ) 


For the subsequent coasting subsegment, the variational equation, with partitioned state 
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transition matrix, takes the same form as in the impulsive formulation, 4 


Sr k - v k 5t k 


A k ,r B kT 


Sr T — v^5t T 

_ -a~5t k _ 


C k ,T D k: T 




( 3 ) 


For this formulation, both the initial and final times of the arc (t k _i and t k ) are fixed, 
though that is not a requirement. The initial position r^_i, velocity v&_i, and mass m k _i 
are also fixed. The mass flow rate, m Sfc _ 1 , is also fixed. It is important to note that vj = 
(and therefore <5vJ = Jv^). Furthermore, — a + (a^ — a j,)StT- 

Incorporating these substitutions, the first two vector variational equations from Equation 
(1) and Equation (3) can be combined to give an expression for dr/-: 


5r k = 



Su t - 1 

(A kjT GT, k -i + B kT J Tk _ i) B k)T { a T — aj) 





5t T 


As in the impulsive Level I method, 4 a minimum norm solution is selected to obtain the 
desired change in the control variables, 


i 

5t T 


M T (MM T )~ 1 6r k , 


(5) 


where 


M = 


(Ak,TGT,k-i + ~ a r) 


( 6 ) 


A minimum norm solution identihes the smallest change in the control parameters, in this 
case <5u^_ 1 and Str, that lead to the desired changes in the constraint errors. Of course, these 
corrections are linear in nature and, as such, an iterative process is required to converge on 
the specified constraints in the nonlinear system. 

To determine an initial guess for the finite burn (i.e. thrust direction, thrust magnitude, 
and burn time), first the impulsive Level I process is used to compute an impulsive correction. 
The impulsive Av direction is used as an initial guess for the thrust direction. The desired 
thrust value (i.e. the thrust of the engine) is used as an initial guess for the thrust magnitude. 
Finally, the burn time is deduced by starting with the rocket equation, 


Av fc 


Isp9 o^(l 


At burn ' 

m k 


(7) 
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substituting in m gk = Uk/I sp g 0 , where u k is the thrust magnitude, and rearranging to obtain 


It is important to note that because this initial guess is based on an impulsive maneuver, 
the terminal error after the first iteration can be very large when the burn duration is long. 
The burn direction is assumed to be constant throughout the entire maneuver, and so small 
errors in direction can be greatly magnified by the end of a long burn. 

1. Controlling Thrust Magnitude 

The 8 u+ 4 in Equation 5 implies changes in both the direction and magnitude of the thrust 
vector, Ufc. For a constant thrust engine, however, a change in the thrust magnitude is 
clearly not desirable. Fortunately, there are several workarounds for controlling the thrust 
magnitude to the desired value. 

The first approach is to simply ignore the change in thrust magnitude suggested by 5u^_ 1 
and only use the change in thrust direction. The new thrust direction along with desired 
thrust magnitude are used to create a new thrust vector in the equations of motion. The 
Level I process converges using this technique, but convergence is slower since the updates 
to thrust vector in the equations of motion do not match the updates suggested by the 
minimum norm solution. 

The second approach for controlling the thrust magnitude is to use a thrust biasing 
technique. This technique is similar to final position biasing commonly used in perturbed 
Lambert targeting. insert FDO Console Handbook citation. The technique uses the same 
Level I process, but implements the full thrust vector update (direction and magnitude) 
suggested by the minimum norm solution in the equations of motion. a At convergence, 
8 r k = 0 , but the converged thrust magnitude will different from the desired value. Let 
A Uk = Udesired ~ u converged ■ The Level I process is then repeated with the same initial 
guess for thrust direction, but the initial guess for the thrust magnitude is biased such that 
^initial (new) = Uinitiai (old) + Aw fc . The initial guess for the burn time is also updated using 
Equation 8 to reflect this change in thrust magnitude. At convergence, the thrust magnitude 
will again be different from, but now much closer to the desired value. This process of biasing 
and re-converging is repeated (typically 3 to 5 iterations) until A u k = 0. 

Figure 4 provides example output from a Level I process demonstrating the thrust bi- 

a Since rh gk is assumed constant, in effect this creates a fictitious variable I sp engine since m gk = u k / T sp go 
but u k is changing after each iteration. This is only temporary, however, since the thrust biasing technique 
ultimately brings the thrust magnitude back to the desired value, and hence I sp ultimately returns to its 
assumed value as well. 



( 8 ) 
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asing technique. This particular example is targeting the Trans-Earth Injection 1 (TEI-1) 
burn as part of the Earth entry targeting process. In Figure 4, the first four lines show the 
convergence of the impulsive corrector to provide the initial guess for the thrust direction. 
The next six lines demonstrate the first pass through the Level I process. Note how the ini- 
tial guess of the thrust magnitude (UMAG) is the desired value of 33361 N. At convergence, 
however, the thrust magnitude is 32802 N - 559 N lower than the desired value. For the 
second pass, the initial guess of the thrust magnitude is biased upward to 33920 N, and at 
convergence results in a value of 33355 N - now only 6 N lower than the desired value. The 
third and final pass through the Level I process converges back to the desired value of 33361 
N. 
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Figure 4. Thrust Biasing in the Level I Process 


A third approach is to add the thrust magnitude as a constraint in Equation 4. Let 


Ctfc U k 


C'l;. ' U;, 


( 9 ) 


Then 


Sa k = u k - u k * 


( 10 ) 
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and 


( 11 ) 


da k _ A T da k _ 
dt T ~ 

Equation 4 then becomes 


1 — 
Sh 


{Ak,TGT,k- 1 + B k pJT,k- 1) 

B kj T (&T ^t) 
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+1 

°o 
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Q 


b Jfi 

<3 
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5tT 


During implementation, however, this method lead to divergence in 8r k . The minimum norm 
solution in this case corrected the thrust magnitude (no matter how large the error) fully in 
one iteration to the desired value and did not deviate in subsequent iterations. This behavior 
was not surprising given that the partial of the constraint with respect to the thrust vector 
is the thrust unit vector itself. Once the thrust magnitude became “fixed,” it may not have 
been possible for the targeter to correct final position discontinuities with changes in the 
burn time and/or direction alone. 

A fourth and final approach for controlling the thrust magnitude is a hybrid of the second 
and third methods. From the second method, the Level I process is run once to convergence 
using the full thrust vector update from the minimum norm solution. A thrust bias is 
determined and the initial guess for the thrust magnitude is updated. At this point however, 
the thrust magnitude constraint from the third method is imposed prior to repeating the 
Level I process. Running Level I once prior to imposing the constraint improves the initial 
guess and should aid convergence. This approach has not yet been implemented and is left 
for future work. If successful, it would require only 2 iterations of the Level I process rather 
than 3 (which is typically the best case scenario) for the thrust biasing technique. 


2. Variable Scaling 


A well-known tool for aiding (and sometimes enabling) the convergence of iterative processes 
such as non-linear targeting is the scaling of both the control variables and constraints to 
0(1). This enables the targeter to evenly adjust control variables and meet constraints of 
varying orders of magnitude. In a physical sense, this is analogous to a person artificially 
adjusting masses and sizes to make juggling a bowling ball, a peanut, and an inflated ballon 
feel more like juggling a tennis ball, a baseball, and a racquet ball. 

To implement scaling in the Level I Process, first note that Equation 4 relates variations in 
the constraints (<5r fc ) to variations in the controls (5u/y, and Stx) through partial derivatives 
as 


dr k dv k 


1 

+1 

_ du +_ x dt T 


5tx 


( 13 ) 
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Prior to forming the minimum norm solution, the controls, constraints, and partial deriva- 
tives are scaled by suitable factors so that 


U fc-l sc z " 

_ u i~ i 

U'scl 

N 

ii 

tx 

tsci 

Sr k 3cl = 

r s d 

dr k 3cl 

dr k 

Uscl 

dr k scl 

_ dr k 

tsd 

du t-l scl 

du k- 1 

7* scl 


dt T 

Tsci 


and therefore, 


dr kscl dr kscl ' 


6u t- isa 

. du t- i scl dt T sc i . 


St T sc i 


The scaled controls are then updated using the minimum norm solution. To begin the next 
iteration, the controls are un-scaled and then used in the equations of motion to determine 
the new constraint variations and partial derivatives. As before, these are re-scaled, a new 
minimum norm solution is found, and the process is repeated until convergence. 

Since the constraints are always being driven to zero in Level I, it is sufficient to set 
r sc i = 1. u sc i can nominally be set to the desired thrust value. t sc i can be set to the 
initial estimate for Atb urn from Equation 8 if the thrust is sufficiently high. If the thrust is 
sufficiently low such that the impulsive approximation is not valid (at least for the initial 
guess), then the burn duration will be much longer and t sc i should be increased appropriately. 

Finally, note that scaling can and should also be implemented in the Level II Process. In 
Level II, for the velocity continuity constraint, the controls are r fc _!, t k _ i, r k , t k , r fe+1 , and 
t k +i, the constraint is Av^, and the partial derivatives are given in Equation 18. Note that 
in Level I, r k was a constraint, and therefore it was sufficient to set r sc i = 1. In Level II, 
however, r fc _ 1; r fc , and r fc+1 are now controls, and must be scaled according to their expected 
magnitudes at convergence. Similarly, in Level I tx was a burn time, and therefore t sc i was 
set to the initial estimate for A tb U m- In Level II, however, t k - i, t k , and t k+ 1 are all patch 
point epochs, and should be scaled according to their expected values at convergence. This 
same approach is also taken for the maneuver sum constraint. 

3. Level 1 Earth Entry Targeting 

The following example shows a comparison of solutions obtained from a Level I process 
based on the impulsive and finite burn targeters. In this simulation, the Level 1 algorithm 
is implemented for a trajectory segment from a given TEI-3 position to Earth entry. In 
keeping with the overall structure of the two-level corrector, only the terminal position 
vector is targeted for this example. TEI-3 occurs on April 6, 2024 10:07:23.8635 UTC, and 
the desired entry state is as follows: 
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• Geodetic Altitude (km) 121.92 

• Longitude (deg) 175.6365 

• Geocentric Azimuth (deg) 49.3291 

• Geocentric Flight Path Angle (deg) -5.86 

Because altitude is a function of position only, the error in meeting the desired altitude re- 
lates directly to the ability of the Level I process to meet the desired final position. Likewise, 
with the final time fixed, the final longitude error can also be judged by the final position 
error. Because the main engines are employed for this scenario, the impulsive assumption 
is reasonably accurate and the results for the finite burn algorithm should therefore match 
those of the impulsive fairly closely. Table 1 lists the convergence data for this case. 


Table 1. Level I Targeting Results 



Impulsive 

Finite Burn 

Initial Position Error (km) 

3.48696e5 

3.91018e6 

Final Position Error (km) 

7.29519e-6 

8.83104e-5 

FPA Error (deg) 

0.0112 

0.0123 

Azimuth Error (deg) 

0.0661 

0.0782 


The maneuver Av and final constraint error values are very similar for the impulsive 
and finite burn Level I processes. It is also evident, from Table 1, that there appears to 
be a significant discrepancy in the initial position error provided to each targeter. The 
discrepancies between the startup arcs for each targeter originate from the estimation of 
the finite burn maneuver previously described. That is, the direction of the burn, which 
is assumed to be inertially fixed, is initially aligned with the Av vector computed by the 
impulsive targeter. For a burn of finite duration, the direction of the burn becomes more 
significant as the integration time increases. Thus, an error in the thrust direction is the 
source of the increased initial position error reported in Table 1. As a result of this larger 
initial error, the finite burn targeter requires added iterations to converge on a solution in 
this case. However, the larger error is also a good indication of the accuracy of the impulsive 
assumption in this example. 

In general, both cases considered here required an increased number of iterations relative 
to later examples. This is attributed to the sensitivities traditionally associated with a 
Level I process. Specifically, the success of a Level I process is sensitive to the integration 
time. Since this particular example employs a Level I process to transfer the vehicle from 
the vicinity of the moon, at TEI-3, to the entry interface at Earth, the time of flight is 
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too long for the number of control parameters available. Here, a Level 11 process 4 becomes 
useful because it allows for an increased number of control parameters and also the ability 
to incorporate an arbitrary number of constraints. The development of a Level 11 process 
that accommodates finite burn maneuvers is discussed next. 

B. Level II Process 

In the classical two-level corrector, 4 velocity discontinuities between coast segments arise due 
to the Level I process. This is also applicable to the finite burn formulation, except at the 
point where a finite burn maneuver is initiated. Here, the burn segment is always assumed 
to start with the same initial velocity as the terminal velocity of the preceding arc. The 
duration of the burn does not exceed the time associated with the following patch point. 
Thus, a velocity discontinuity can occur, during the Level I process, at the point where the 
coast subarc, as defined in Figure 3, joins with the following trajectory segment. Although 
this problem at first seems identical to the impulsive maneuver targeting, since the velocity 
discontinuity falls between two coast arcs, the partial derivatives for dvjT with respect to 
5r k _i, t k _ i, 5r k , and t k differ due to the thrust segment at the beginning of the arc. 

Recall from the Level I formulation that = vj at the terminal point of the burn 
arc and thus that <5v^ — — a ^5tr + (aj^ — a ^)5tT- For the Level II process, 

b'm k _i = rhg k _ 1 = 0 and u^_, = 0. It is still assumed that fn g is a fixed constant, i.e. 

= 0. Using these relationships and assumptions, along with Equations (1) and (3), 
an expression is found for in terms of the state at patch point k — 1 and the state 
transition matrix, 

= Ck,T[AT,k-l(fi r k-l — 

+ H Tjfc _ 1 (dvj|'_ 1 - a^_ 1 5t fc _ 1 ) + + G T ,k- 

(15) 

+ D ki T[CT,k-l(fck-l ~ + DT,k-l(^ k -i ~ 

+ + J T ,k- 1^_ 1 + (a^ - aJ)df T ] + a ^5t k . 

In order to write only in terms of the Level II control variables, the first vector equation 
from Equation (1) is used to solve for 5v^_ 1 , and 5tx in terms of those control 
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variables. From the minimum norm solution, 


5v fc-i 

Su t - 1 

St T 


Z T (ZZ T ) 1 [5r k -v k 5t k 


— + B kt TCT,k-l)(^k-l ~ V^'_ 1 5tfc_ i) 
+ {A kT B T ^k-l + 

— (Ak,TET,k-i + ^fc,r-f^T,fc-i)hr^ c _ 1 htfc_i]. 


( 16 ) 


where Z = 


(^ 4 fc,T^T,fc-i + Bk^Dr^k- 1) (Ak,TGT,k-i + B k ^JT,k- 1) B k ,T( &r ~~ a r) 


With 


this expression, the partial derivatives of Av^ with respect to each control variable can be 
found using the same method as in the impulsive formulation. Let 


Z = 


(C^tBt ^- i + D kT D Tk _ i) (C kj TGT t k - i + D kT J Tj k-i) D kT (c 


— a. 


Z T (ZZ 


T \- 1 


(17) 

Then, because it is assumed that the arc from patch points k to k + 1 is a coast arc, the 
partial derivatives of Av k are 


<9A v fc 
<9r fc _i 
dAv k 
dt k - 1 


dAv fc 

<9r fc 
dAv k 
dt k 
dAv k 
dr k+ i 
dAw k 
dt k -\~ i 


— [(C’fc.r^T.fc-i + D k ^CT,k- i) — Z(Ak,TAr,k - i + 

— ([(Cfc,T^T,fc-l + D kj THT,k- 1) _ Z(A k ^ET,k-l + -Sfc,T^T,fc-l)]hT.^__ l 

— [(Cfc,T-BT,fc— i + Dk^D^k- i) — Z(A kt TBT,k-i + 5fc i T£ ) T,fc-i)]a^"_ 1 

— [(Cfc.T^T.fc-i + D k^pCr ,k—i) — Z(A k ^AT : k-i + Bk^Cr^- i)]v^_ 1 ), 

S fci,jA+i,* v fc + a+ - (a" - Zv“), 

D — 1 

^£:+l,AO 
— -®fc+i,fc v fc+r 


(18) 


The above partials are employed in the standard level If process. 4 


C. Maneuver Sum Constraint 

In addition to the velocity continuity constraint, endpoint and interior path constraints may 
be imposed during the Level If process. 16 One such constraint is on the total Av sum of 
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the maneuvers. The finite burn formulation of this constraint is based on the impulsive 
maneuver sum constraint. 17 Only the composition of the associated partial derivatives and 
the error calculation changes. 

To derive the burn maneuver constraint, it is necessary to determine the partial deriva- 
tives of the magnitude of Av*,, i.e. the maneuver that results from the burn at patch point 
k, with respect to the Level If control variables. From the rocket equation, Av^. is given by 

Avj; = -I, r g 0 ln( 1 - (19) 

m k 


where A t burn = 
k is given by 


t T — t k . The partial derivative of Av fc with respect to A at patch point 


dAv k 
dAt 


burn 


Isp9o( 


m k 


)( 


m 


9k 


Wlk TTl g k /\.tb urn Tflk 


( 20 ) 


Next, the partial derivatives of A tb urn with respect to the control variables are necessary. 
This is determined using the variational equations from points k — 1 to k, k to T (the 
termination of the burn segment), and k + 1 to T. Recalling that A tburn — t t — t k , the 
partials are found to be 


dAt burn 
<9r fc _i 

9At burn 

dt k ~ i 

9 At burn 
dr k 

dAtburn 

dt k 


d Atburn 

dr k+1 
9 At f )Ur n 


U k 

|Aa 7 


ut 


I Aa 7 


ut 


||Aa r || 


^D T , k B~\ k - SB T , k B~\ k \ 
(D T)k B~\ k - S B Tjk B k \ k )v k _ v 


\{CT, k + DT, k D k ^ k _iB kk _k) — S(AT )k + BT, k D ktk -iB kk _ 1 )], 


ut 


I Aa 7 


([(C’t.a: + DT,k.D ktk -iB kk _ l ) — S(AT )k + BT, k D ktk _iB kk _k)\ 


( Dr,k — SB T!k )(a k — a k ) — (H Tk — SE Ttk )m g ), 


ut 


I Aa 7 


-(O 


T,k + 1 




T,k+ l)i 


Ut 


dt 


fc+1 


l|Aa r || 


[{Ct, fe+i — SAx }k + i)v fc+1 + (Dx,k+ i — SBx )k + i) a fc+1 ], 


( 21 ) 


where S = — Gx k B Tk+1 and S = —Jxk Dx k +i S' T (S'S' T )~ 1 . These partials leave 

one term, containing 5m k , unaccounted for. The cost of a maneuver at patch point k depends 
on the duration of the burn at k and the initial mass m k . Since m g is a fixed, constant value, 
Sm k depends only on the previous burn durations. Thus, the initial mass at the beginning of 
a maneuver will have a dependence on the positions and times associated with any previous 
maneuvers that have occurred, as shown above. Using the chain rule, the final form of the 
partial derivative of the constraint a (the sum of all the burn Avs) with respect to any 
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control variable (3 k in the set of control variables associated with patch point k is 


da 

d/3 k 


JV Av 

X 


n = 1 


OAv n 

d/3 k 


( 22 ) 


where N& v is the total number of maneuvers implemented along the trajectory and 

dAv n _ dAv n dA t burnn dAv n dm n 

df3 k dAt burnn df3 k dm n df3 k ' 


Because the mass at the time of a burn, m n , depends on the propellant mass consumed 
during the previous burns, 


dm n diri n —i . dAt burr ^ ^ 

— — 777 

Sft Sft 9/ i 


(24) 


A similar relationship exists for the remaining mass partials (m n _ 2 to mi) with respect to 
/3i. These partials are then employed during the Level 11 process. 4 


IV. Simulation and Results 

The two level targeting algorithm is applied in the following section in both the impulsive 
and finite burn configurations. This is done in order to demonstrate its performance under 
both configurations and to gain insight as to applications of each. Then, as a final metric, 
an optimal trajectory is generated to use for performance comparisons. 

For each case, the same initial conditions of the lunar orbit will be used as follows: 

• Epoch: 4-Apr-2024 15:30:00 TDT 

• Initial mass: 20339.9 kg (total fuel = 8063.65 kg) 

• Main Engine Thrust: 33,361.6621 N 

• Main Engine Isp: 326 sec 

• Auxiliary Engine Thrust: 4,448.0 N 

• Auxiliary Engine Isp: 309 sec 

• State (J2000 Moon-centered inertial frame): 

- X: -1236.7970783385588 km 

- Y: 1268.1142350088496 km 
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- Z: 468.38317094160635 km 

- Vx: 0.0329108058365355 km/sec 

- Vy: 0.589269803607714 km/sec 

- Vz -1.528058717568413 km/sec 

Likewise, the same terminal target conditions will be used for each case as shown below: 

• Geodetic Altitude (km) 121.92 

• Longitude (deg) 175.6365 

• Geocentric Azimuth (deg) 49.3291 

• Geocentric Flight Path Angle (deg) -5.86 

A. Finite Burn Example with Main Engines 

The first case is representative of a nominal Earth return during which the maneuvers are 
performed by the CM main engine. The initial guess hie consists of 12 ’’patch points” or 
states taken from an optimized trajectory. The first patch point corresponds to a state on 
the initial lunar parking orbit. The interior patch points correspond to the states and epochs 
at each of maneuver locations (TEI-1, 2, 3 and TCM 1, 2, 3) and some additional waypoints 
along the trajectory. The final patch point in the initial guess is the state and epoch at the 
desired entry interface (El). 

For this case, both the impulsive targeter and the finite burn targeter are executed in 
order to find a feasible trajectory that satisfies the specified terminal constraints, while 
keeping the Av sum of the individual maneuvers within the available fuel budget. Table 
2 compares the individual maneuvers and final Av sum for the impulsive solution and the 
finite burn solution. The burn parameters for each finite burn maneuver are given in Table 
3, and final constraint errors are listed in Table 4. 


Table 2. Maneuver Data 


Maneuver 

Impulsive Av (km/s) 

Finite Burn Av (km/s) 

TEI-1 

0.6619 

0.7348 

TEI-2 

0.3257 

0.2561 

TEI-3 

0.4115 

0.4087 

Total 

1.3991 

1.3996 


The individual maneuvers and total Av sum for the finite burn targeter are fairly similar 
to the impulsive targeting results, which is to be expected given that the burn durations 
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Figure 5. Initial Guess and Impulsive Solution 


Table 3. Burn Data 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

TEI-1 

400.2308 

4175.128 

TEI-2 

119.2402 

1243.890 

TEI-3 

171.6074 

1790.174 

Total 

691.0784 

7209.192 


Table 4. Constraint Error Data 



Impulsive Algorithm 

Finite Burn Algorithm 

Iterations 

20 

6 

Altitude (km) 

-1.0e-8 

-5.2e-8 

Flight Path Angle (deg) 

2.8e-10 

-3.3e-10 

Longitude (deg) 

4.0e-8 

-5.0e-7 

Flight Path Azimuth (deg) 

-2.0e-7 

1.6e-7 
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with the main engines are short enough for an impulsive assumption to be used. It should 
be noted, though, that for this particular case, the impulsive algorithm requires several more 
iterations to converge than the finite burn algorithm does. This suggests that the impulsive 
assumption, while still valid, may be reaching its limit. 

B. Finite Burn Example with Auxiliary Engines 

For this example, a main engine failure is assumed to occur after TEI-1 and the auxiliary 
engines are used to perform the final two maneuvers. Figures 6 through 8 show the closeup 
views of TEI-1, TEI-2, and TEI-3, and the solid red portions of these closeups indicate the 
segments of the trajectory in which the engines are thrusting. Burn data for each maneuver 
is listed in Table 5, and the final constraint errors are given in Table 6. 



Figure 6. Finite Burn Solution with Auxiliary Engines, TEI-1 Closeup 


Table 5. Burn Data Using Auxiliary Engines 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

Av (km/ s) 

TEI-1 

363.5548 

3792.531 

0.6255 

TEI-2 

949.2614 

1392.946 

0.2666 

TEI-3 

1400.5756 

2055.205 

0.4418 

Total 

2713.3918 

7240.682 

1.3339 


The finite burn algorithm is able to meet the entry and cost constraints in only a few 
more iterations than the impulsive targeter required. Interestingly, the total Av for this case 
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Figure 7. Finite Burn Solution with Auxiliary Engines, TEI-2 Closeup 



Figure 8. Finite Burn Solution with Auxiliary Engines, TEI-3 Closeup 
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Table 6. Constraint Error Data with Auxiliary Engines 


Iterations 

23 

Altitude (km) 

4.1e-7 

Flight Path Angle (deg) 

4.5e-10 

Longitude (deg) 

2.1e-6 

Flight Path Azimuth (deg) 

-2.6e-7 


is about 0.065 km/s lower than the impulsive solution. This will not always be the case; 
numerous feasible solutions can exist for any given set of patch points. For this particular 
example, the finite burn Level II process identified a lower cost solution than that determined 
with the impulsive targeter. In both cases, the total cost constraint is always enforced to 
ensure that the total cost is within the available fuel budget. 

C. Optimized Finite Burn Trajectory 

As a final step in this analysis, the trajectory generated from the finite burn targeter with 
auxiliary thrusters was optimized using Copernicus. The results of this optimized run are 
available in Table 7. The total Av for the optimal run is 1.2413 km/s. This is an improvement 
of approximately 0.09 km/s of Av over the finite burn targeting solution. However, the cost 
constraint imposed during the Level II process specified the total cost should not exceed 1.40 
km/s. Specifying a lower boundary on this constraint may have identified a similar solution. 
It is always important to bear in mind that a targeter does not seek optimal solutions, only 
feasible solutions. If a feasible solution exists in the vicinity of the initial guess, either the 
impulsive or the finite burn targeting algorithm, can typically identify it. 


Table 7. Optimal Burn Data Using Auxiliary Engines 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

Av (km/s) 

TEI-1 

3327.008 

3501.37 

0.6040 

TEI-2 

873.040 

1504.85 

0.2993 

TEI-3 

1312.285 

1828.22 

0.4059 


D. Finite Burn Example Over The Lunar Cycle 

To further test the finite burn algorithm, return trajectories were generated over several days 
spanning the lunar cycle from February 1-28, 2024, 0:00:00 TDT, again using the auxiliary 
engines for the second and third TEI maneuvers. For these runs, only two entry constraints 
are targeted: 
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• Geodetic Altitude (km) 121.92 

• Geocentric Flight Path Angle (deg) -5.86 

In this simulation, the input patch points to the finite burn algorithm come from a 
converged impulsive trajectory with the same initial point and entry targets. The total Av 
of the impulsive trajectory for each case is 1.50 km/s. Table 8 lists results for days 1, 3, 6, 
10, 13, 16, 19, 22, 25, and 28 of the lunar cycle. Figures 9(a)-9(c) show some examples of 
these trajectories. 


xio* xio* 



Figure 9. Return Trajectories (MCI) at Day 1 (a), Day 13 (b), and Day 28 (c) of the Lunar 
Cycle 

With the exception of the simulations for February 6 and February 16, the targeter 
converges to valid solution within the entry and cost constraints for each day tested over the 
lunar cycle. The solutions for February 6 and 16 satisfy the entry constraints, but would 
not converge when the total cost constraint of 1.50 km/s was imposed. The values listed in 
Table 8 are the total cost of the converged solution without the maneuver sum constraint 
active. The precise effect of the quality of initial guess on a given date will be the subject of 
future investigations. 
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Table 8. Burn Data over the Lunar Cycle 


Day 

TEI-1 

TEI-2 

TEI-3 

Total Cost 


AV (km/s) 

Duration (s) 

AV 

Duration 

AV 

Duration 

Total AV 

Iterations 

1 

0.5963 

348.1673 

0.4741 

1648.5129 

0.3995 

1202.2754 

1.4698 

4 

3 

0.6837 

393.6876 

0.4393 

1492.3184 

0.3585 

1067.4428 

1.4814 

4 

6 

0.5901 

344.8973 

0.5271 

1820.9042 

0.3929 

1165.7458 

1.5101 

15 

10 

0.5849 

342.1671 

0.5808 

1992.6327 

0.3259 

961.9933 

1.4916 

4 

13 

0.5869 

343.1966 

0.3661 

1299.5151 

0.5079 

1561.4187 

1.4609 

4 

16 

0.6806 

392.0921 

0.3205 

1111.0622 

0.6251 

1856.5447 

1.6262 

39 

19 

0.5778 

338.3921 

0.5925 

2033.8077 

0.3189 

940.8576 

1.4892 

4 

22 

0.7004 

402.2770 

0.3915 

1332.7907 

0.4080 

1217.7097 

1.4999 

11 

25 

0.5862 

342.8320 

0.6318 

2149.2973 

0.2740 

801.4648 

1.4919 

4 

28 

0.5861 

342.7945 

0.4912 

1709.0768 

0.3882 

1167.9040 

1.4656 

4 


E. Delayed Patch Point Simulations 

Another test of the finite burn algorithm is whether or not it can converge on a feasible 
solution given a set of patch points that are not current. A set of patch points corresponding 
to a current or future departure time may not always be available, especially when ground 
communications are lost. The algorithm must therefore be able to converge even when the 
departure time listed in the input hie has already passed. For this example, the input patch 
point hie from the February 1 run in the previous section is used. As before, the auxiliary 
engines perform the TEI-2 and TEI-3 maneuvers. To ensure that the characteristics of the 
initial lunar orbit remain the same, the patch points are converted to the MCI frame before 
the time delay is introduced. The initial time of the simulation is perturbed hrst for only 3 
hours, then for a full 12 hours. Figures 10 (a)- 10(b) show the difference in the initial and hnal 
converged trajectories for time delays of 3 hours (10(a)) and 12 hours (10(b)). The initial 
trajectory appears as a dashed line in Figures 10(a)-10(b). The initial entry constraint errors 
due to the time delay are given in Table 9, and the final constraint errors for the converged 
trajectories are listed in Table 10. 


Table 9. Initial Entry Constraint Errors 


Delay (hr) 

Altitude Error (km) 

FPA Error (deg) 

3 

1.1040e4 

1.6464 

12 

4.5332e4 

10.4268 


Even with a 12 hour delay, the finite burn algorithm is able to End a feasible solution 
that satisfies the entry constraints and fuel budget. The 12 hour delay trajectory even shows 
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(a) 3 Hours 


(b) 12 Hours 


Figure 10. Initial and Final Trajectories (MCI) at 3 Hours (a) and 12 Hours (b) Delay 


Table 10. Convergence with Delays 


Delay (hr) 

Total AV (km/s) 

Altitude Error (km) 

FPA Error (deg) 

Iterations 

3 

1.4655 

-2.8218e-5 

-9.4476e-8 

10 

12 

1.2961 

6.0979e-6 

1.2184e-8 

12 


a surprisingly marked decrease in the total Av cost for this case. As with the low-cost 
auxiliary engine solution discussed previously, it is the sensitive dynamics of the system that 
accounts for this result rather than any special property of the finite burn targeter. These 
results do, however, underscore the flexibility of the two-level targeting structure; instead 
of trying to match a previously determined nominal trajectory, the algorithm explores the 
nearby solution space and is able to converge on a trajectory that is more desirable than the 
one suggested by the initial input. 


V. Conclusions 

This paper presents a two-level targeting algorithm for finite burn maneuvers. The algo- 
rithm is adapted from a classical impulsive two-level targeting algorithm. The development 
of the finite burn versions of both the Level I and Level II process are discussed, and a 
total mission cost constraint, originally developed for the impulsive algorithm, is modified 
and adapted to the finite burn problem. The algorithm keeps the same structure and much 
of the simplicity of the impulsive two-level corrector even though the partial derivatives re- 
quired for the calculations are far more complex. For testing, the algorithm is applied to the 
trans-Earth injection phase of the Orion mission. Results are compared for three different 
cases: using the impulsive algorithm, using the finite burn algorithm with main engines for 
all three maneuvers, and using the finite burn algorithm with only the auxiliary engines 
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for the final two maneuvers. These results show that the finite burn algorithm is able to 
converge on a feasible solution even in the case of a main engine failure following TEI-1. Due 
to the additional complexity of the finite burn model, the algorithm does exhibit increased 
computational overhead in contrast to the impulsive targeter. However, the present formula- 
tion addresses the need for a targeting algorithm that accommodates the main engine failure 
scenario while meeting all the specified constraints. 
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